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ABSTRACT 

Fluid-structure interaction (FSI) is the commerce between a flexible structure and the surrounding fluid 
gives a variety of phenomenon. Application of this phenomenon in many areas like vibrational analysis of aircraft 
wings, machinery, bridges, etc. This computational analysis leads to many aerodynamic and structural dynamic 
problems. In the present article, we review the structural analysis of cropped delta wing aircraft with fence and 
without the fence. One of the key aspects of developing a better design from the data collected by FCI. Several such 
techniques are reviewed in this paper. As a part of this work, the experimental model setup for analysis the cropped 
delta wing by using wind tunnel and validate with software analysis values. The model ofa cropped delta wing with 
the fence at different Mach number and angle of attack are selected to high spot the futuristic of computational 
problems. 
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1. INTRODUCTION 

The wing of an aircraft has a major role in the aircraft, contributing all its lift produced to sustain in the air. 
The entire load of the aircraft including all the passenger, fuel, cargo, Artillery if it is a military aircraft, is been 
carried by the wings. So a tremendous amount of load is acting on the aircraft wing structure at all operational 
levels. To maintain its all-important aerodynamic shape, a wing must be designed and built to hold its shape even 
under extreme stress. Basically, the wing is a framework composed chiefly of spars, ribs, and (possibly) stringers 
[25]. 
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Figure 1: An Internal Structure of an Aircraft Wing 


Spars are extended lengthwise of the wing (crosswise of the fuselage). In fiight, the force of the air acts against 
the skin. From the skin, this force is transmitted to the ribs and then to the spars. The ribs are the parts of a wing which 
support the covering and provide the airfoil shape. Their purposes are providing shape, bearing fiight stress, and these are 
called compression ribs [18]. 

It is a comparative study between the deformations due to the loads over the modified wing with the original. So 
the loads due to weight remain constant and only the aerodynamic loads have to be taken for the analysis. 
The Aerodynamics loads are obtained from the CFD analysis performed and results can be directly used in the structural 
analysis [4]. The most modern Computer Aided Engineering (CAE) technique called Fluid-Structure Interaction (FSI)[5] is 
applied. It is the interaction of some movable or deformable structure with an internal or surrounding fluid flow. 

The second phase of this work is the testing of prototype model by using the wind tunnel. During a test, the model 
is placed in the test section of the tunnel and air is made to flow past the model. In some wind tunnel tests, flow 
visualization techniques are used to provide diagnostic information about the flow around the model [23]. Smoke is used to 
visualizing the flow that is away from the surface of the model. Smoke can be used to detect vortices and regions of 
separated flow. It can be injected from the surface or dispersed with a hollow wand that can be moved through the flow 
field. The disadvantage of smoke is that it does not work well at higher speeds (greater than -300 mph), the smoke must be 
introduced at the proper location without altering the flow, and the smoke can leave a residue in the tunnel or on the model, 
depending on the type of smoke employed. 

2. MATERIAL SELECTION 

The main doubt about the metallic wing was regarding the kind of aluminum to be used in its components. In the 
metallic projects of conventional aviation there is a concern with the final price of the aircraft, therefore, in the areas of the 
components with less structural load, the aluminum used is of a less noble alloy and, thus cheaper. The mechanical 
properties of the aluminum depend on the width of the plate and orientation of the grains. For this work, in order to 
maintain the use of only one material, even when there is the reduction of the width of the skin or of the stringers 
throughout the span, it was decided to adopt a plate with an average width of 0.25” (6.35mm) to the whole model [4]. The 
orientation of the grains on the plate also infiuence the properties, the longitudinal orientation of the grains was chosen for 
it demonstrates a greater endurance than the others. The elastic isotropic properties and the mechanical properties of the 
Aluminum used are exposed in Table 1. 
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Table 1: Mechanical Properties of Aluminum 


E 

(MPa) 

Gi2 

(MPa) 

Ni 12 

F t x 

(MPa) 

Ftv 

(MPa) 

Fcx 

(MPa) 

F C y 

(MPa) 

701016 

26889 

0.31 

489 

413 

420 

282 


3. DESIGN AND FLUID STRUCTURE INTERACTION (FSI) ANALYSIS 

The dimensions of the structure consist of a root chord of 8474.87 mm, a semi-span of 5500mm and a tip chord of 
550.87 mm (Figure 2). 



Figure 2: Spars and Ribs of a Wing 

The Delta wing structure created using CATIA V5R21, was subjected to aerodynamic loads that are obtained 
from the CFD analysis [16]. The entire wing structure was constructed using 3mm sheet metal. The entire component 
created without any trimming in order to reduce the complexity of meshing and analysis. The four spars starting from the 
root rib and converging towards the tip rib. 

The 12 ribs were placed uniformly over the span leaving a distance of 500mm between two ribs. The first spar 
was placed at 2000mm from the leading edge of the root rib. And next two spars placed at 2000mm from the nearest spar. 
The fourth and last spar placed at a distance of 2200mm from the third spar, equally carrying all the loads. A sheet metal 
skin of 3mm thickness was given to the entire structure. 



Figure 3: Wing Structure in Catia Software 


3.1. Meshing of the Wing Structure 

The meshing was done in three steps. In the first step, the material Aluminum alloy was applied from the ANSYS 
Workbench Engineering Data section. Then the Structure was assembled into a single unit. In the final step, the 
Mechanical (Model) option in ANSYS Workbench was used to create a finite model. The maximum element size was 
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given as lOOmm, the number of nodes are 287983 and number of elements are 140135. 



1 .500 


Figure 4: Wing Structure Meshing in ANSYS 

Here the structure and CFD models are created separately and solved using Multi-physics. The results from CFD 
analysis are transferred into the structure analysis. The aerodynamic loads are imported as pressure and solved to get 
various results like total deformation, equivalent Stress, equivalent Elastic Strain and directional Deformation [24]. 

The FSI analysis was carried out at 0°, 9°, 18° angles of attacks. The details of finite element model are given in 

Table 2. 


Table 2: Finite Model Description 


Object Name 

Part 1 

Part 2 

State 

Meshed 

Meshed 

Stiffness Behavior 

Flexible 

Flexible 

Coordinate System 

Default Coordinate System 

Default Coordinate System 

Reference Temperature 

By Environment 

By Environment 

Material 

Assignment 

Aluminum Alloy 

Aluminum Alloy 

Nonlinear Effects 

Yes 

Yes 

Thermal Strain Effects 

Yes 

Yes 

Bounding Box 

Length X 

8.4765 m 


Length Y 

0.85731 m 


Length Z 

5.5 m 

5.5033 m 

Properties 

Volume 

7.7715e-002 m 3 

4.6761e-002 m 3 

Mass 

215.27 kg 

129.53 kg 

Centroid X 

5.5737 m 

4.618 m 

Centroid Y 

6.3928e-002 m 

9.0979e-002 m 

Centroid Z 

1.8902 m 

1.4773 m 

Moment of Inertia Ipl 

200.2 kg-m 2 

124.34 kg-m 2 

Moment of Inertia Ip2 

894.56 kg-m 2 

736.35 kg-m 2 

Moment of Inertia Ip3 

707.32 kg-m 2 

620.04 kg-m 2 

Statistics 

Nodes 

206825 

81158 

Elements 

103496 

36639 


3. 2. Structural Analysis of Wing without Fence 

The figure 5 shows the loads imported from the CFD analysis. The minimum imported pressure was 0.3483 Pa 
that acts along the leading edge. The maximum pressure was 2225.5 Pa acting at the tip and at the wing fuselage 
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attachment part. 



1 .500 


Figure 5: Pressure Distribution over the Wing Structure 


The figure 6 shows the Total deformation of the wing structure under loading. The unit given is in the meter. 
The minimum deformation at the root given as 8.37xl0‘ 4 m. Along the lateral axis, the deformation increase and the 
maximum deformation occurs at the tip. The maximum deformation is noted to be 7.539xl0‘ 3 m. 



Figure 6: Total Deformation Over the Wing Structure 


The figure 7 shows the Directional deformation of the wing structure under loading. The unit given is in the 
meter. The minimum deformation at the root given as -8.564xl0' 5 m. Along the lateral axis, the deformation increase and 
the maximum deformation occurs at the tip. The maximum deformation is noted to be 2.558xl0‘ 4 m. 



Figure 7: Directional Deformation over the Wing Structure 


The figure 8 shows the Shear Elastic Strain of the wing structure under loading. The minimum deformation at the 
root given as -2.7449x10" 4 . Along the lateral axis, the deformation increase and the maximum deformation occurs. 
The maximum deformation is noted to be 1.8921xlO' 3 . 
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B: Static Structural 

Shear Elastic Strain 

Type: Shear Elastic StrainCXY Plane) 
Unit: m/m 

Olobal Ooordinate System 
Time: 1 

4/8/201 3 11 :57 AM 


0.0018921 Max 

0.001 651 3 
0.001 41 06 
0.001 1 699 
0.0009291- 
0.0006884 
0.00044769 
0.00020696 
-3.3765e-5 s 

-0.00027449 Min 



J2* 


Figure 8: Shear Elastic Strain of the Wing 


3. 3. Structural Analysis of Wing with Fence 


The figure 9 shows the loads imported from the CFD analysis. The minimum imported pressure was 0.24186 Pa 
that acts along the leading edge. The maximum pressure was 2431.02 Pa acting at the tip and at the wing fuselage 
attachment part. 



1 .500 


Figure 9: Pressure Distribution over the Wing Structure 


The figure 10 shows the Total deformation of the wing structure under loading. The unit given is in the meter. The 
minimum deformation at the root given as l.lOxlO’ 3 m. Along the lateral axis, the deformation increase and the maximum 
deformation occurs at the tip. The maximum deformation is noted to be 9.944xl0' 3 m. 



Figure 10: Total Deformation over the Wing Structure 


The figure 11 shows the Directional deformation of the wing structure under loading. The unit given is in the 
meter. The minimum deformation at the root given as -9.755 lxl0' 5 m. Along the lateral axis, the deformation increase and 
the maximum deformation occurs at the tip. The maximum deformation is noted to be 3.2203xl0' 4 m. 
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B: Static Structural 

Directional Deformation 

Type: Directional Deformation(XAxis) 

Unit: m 

Global Coordinate System 
Time: 1 

4/8/201 3 2:05 PM 


0.00032203 Max 

0.00027541 
0 00022879 
0.00018217 
—| 0.00013555 
8.8928e-5 
4.2308e-5 
-4.311 4e-6 
-5.0931 e-5 
-9.755 1e-5 Min 



Figure 11: Directional Deformation over the Wing Structure 


The figure 12 shows the Shear Elastic Strain of the wing structure under loading. The minimum deformation at 
the root given as -3.118xl0‘ 4 . Along the lateral axis, the deformation increase and the maximum deformation occurs. The 
maximum deformation is noted to be 2.23lxl0‘ 3 . 



4. STATIC ANALYSIS OF FSI RESULT 


18°. 


The following tabulation depicts the comparative results of FSI analysis at three different angles say 0°,9°, and 


Table 3: Static Analysis of Wing without the Fence for an Angle of Attack 00 


Results 

Object Name 

Total 

Deformation 

Equivalent 

Stress 

Equivalent Elastic 
Strain 

Directional 

Deformation 

Minimum 

0. m 

1.154e-010 Pa 

0. m/m 

-5.3546e-005 m 

Maximum 

4.3058e-003 m 

5.7819e+007 Pa 

8.1435e-004 m/m 

1.5588e-004 m 

Minimum Occurs On 

Part 2 



Part 1 

Maximum Occurs On 

Part 1 

Part 2 


Part 1 


Table 4: Static Analysis of Wing with Fence for an Angle of Attack 00 


Results 

Object Name 

Total 

Deformation 

Equivalent 

Stress 

Equivalent 
Elastic Strain 

Directional 

Deformation 

Minimum 

0. m 

4.444e-011 Pa 

0. m/m 

-6.8395e-005 m 

Maximum 

6.7434e-003 m 

7.434e+007 Pa 

1.047e-003 m/m 

2.2553e-004 m 

Minimum Occurs On 

Part 2 



Part 1 

Maximum Occurs On 

Part 1 

Part 2 


Part 1 
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Table 5: Static Analysis of Wing without Fence for an Angle of Attack 90 


Results 

Object Name 

Total 

Deformation 

Equivalent 

Stress 

Equivalent 
Elastic Strain 

Directional 

Deformation 

Shear Elastic 
Strain 

Minimum 

0. m 

5.3914e-011 Pa 

0. m/m 

-8.564e-005 m 

-2.7449e-004 m/m 

Maximum 

7.5391e-003 m 

8.808le+007 Pa 

1.2406e-003 m/m 

2.558e-004 m 

1.8921e-003 m/m 

Minimum Occurs On 

Part 2 



Part 1 

Part 2 

Maximum Occurs On 

Part 1 

Part 2 


Part 1 

Part 2 


Table 6: Static Analysis of Wing with Fence for an Angle of Attack 90 


Results 

Object Name 

Total 

Deformation 

Equivalent 

Stress 

Equivalent 
Elastic Strain 

Directional 

Deformation 

Minimum 

0. m 

1.1721e-010 Pa 

0. m/m 

-9.755 le-005m 

Maximum 

9.944e-003 m 

1.0443e+008 Pa 

1.4708e-003 m/m 

3.2203e-004 m 

Minimum Occurs On 

Part 2 



Part 1 

Maximum Occurs On 

Part 1 

Part 2 


Part 1 


Table 7: Static Analysis of Wing without Fence for an Angle of Attack 180 


Results 

Object Name 

Total 

Deformation 

Equivalent 

Stress 

Equivalent 
Elastic Strain 

Directional 

Deformation 

Minimum 

0. m 

7.4067e-011 Pa 

0. m/m 

-9.9138e-005 m 

Maximum 

8.8038e-003 m 

9.8526e+007 Pa 

1.3877e-003 m/m 

2.9003e-004 m 

Minimum Occurs On 

Part 2 



Part 1 

Maximum Occurs On 

Part 1 

Part 2 


Part 1 


Table 8: Static Analysis of Wing with Fence for an Angle of Attack 180 


Results 

Object Name 

Total 

Deformation 

Equivalent 

Stress 

Equivalent 
Elastic Strain 

Directional 

Deformation 

Minimum 

0. m 

6.7987e-011 Pa 

0. m/m 

-1.2279e-004m 

Maximum 

1.2733e-002 m 

1.2988e+008 Pa 

1.8292e-003 m/m 

4.0508e-004 m 

Minimum Occurs On 

Part 2 



Part 1 

Maximum Occurs On 

Part 1 

Part 2 


Part 1 


5. EXPERIMENTAL SETUP AND FLOW VISUALIZATION 

For smoke visualization, the prototype was modeled with the following dimensions in CATIA V5R21. 
The CATIA drawing file was converted into step file, for the prototype to be machined. The method used was Synthetic 
Laser Sintering (SLS) process, a well-known rapid prototyping technique. An additive manufacturing layer technology 
SLS involves the use of a high power laser (for example, a carbon dioxide laser) to fuse small particles 
of plastic, metal (direct metal laser sintering), ceramic, or glass powders into a mass that has a desired three-dimensional 
shape. The laser selectively fuses powdered material by scanning cross-sections generated from a 3-D digital description of 
the part (for example from a CAD file or scan data) on the surface of a powder bed. 
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Figure 13: Designed Model of the Prototype in Catia V5R21 


After each cross-section is scanned, the powder bed is lowered by one layer thickness, a new layer of material is 
applied on top, and the process is repeated until the part is completed. Because finished part density depends on peak laser 
power, rather than laser duration, an SLS machine typically uses apulsed laser. The SLS machine preheats the bulk 
powder material in the powder bed somewhat below its melting point, to make it easier for the laser to raise the 
temperature of the selected regions the rest of the way to the melting point. The material used for prototype manufacturing 
is polyamide. Then the model is fixed into the wind tunnel with the following specification. 


Table 9: Wind Tunnel Specifications 


Parameter 

Specifications 

Test Section Size 

600mm x 600mm x 4000mm. 

Maximum Speed 

45m/sec 

Fan 

Axial Flow fan of Diameter: 1.3 meter 
Maximumrpm: 1500 

Number of Blades: 12 

Hub Diameter: 500mm 

Contraction Ratio 

9: 1 

Contraction length 

1.8m 

Settling chamber 

1.8mxl.8m 

Entry section 

Bell mouthed entry 

Honey Comb Size 

5 Ommx 5 0mmx45 Omm. 

Screens 

Two screens 8mesh and 16mesh stainless steel 

Power 

22 KW / 30HP AC motor, with speed control drive 


5.1. Smoke Visualization at Various Angles of Attacks 



Figure 14: Smoke flow Visualization at an Angle of 
Attack90 (without Fences) 
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Figure 15: Smoke flow Visualization at the Angle of 
Attack 180(without Fences) 

During the test, it has been found that the flow separation occurs at the angle of attack of 18° as shown in figure 
15. 

It was observed that the flow separation for delta wing with fences only occurs after the angle of attack of 19.5°. It 
was also found that the delta wing with fences possesses higher stall angle than the delta wing without fences. 



Figure 16: Smoke flow Visualization at the Angle of 
Attack 90 (with Fences) 



Figure 17: Smoke flow Visualization at an Angle of 
Attack 190 (with Fences) 

5. 2. Aerodynamic Forces 

Six Component strain gauge balance has been used to determine the forces such as lift, drag, pitching moment, 
rolling moment, side forces and yawing moments. The apparatus set up is done as per figure 18. 

The setup is made and is placed in the wind tunnel and is tested at the velocity of 30m/sec at various angles of 
attack. The following table shows the lift and drags comparison at various angles of attack with and without wing fences. 
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Figure 18: Experimental Setup for Force Determination 

It can be noted from the graphs that even with a very slight increase in drag, the increase in lift of wing with 
fences is much higher than that without it and has higher stall angle. 


Table 10: Comparison of Forces in with and without the 
Fence for Various Angles of Attack 


Angle of Attack 

Lift 

Dra, 

j 

Without Fence 

With Fence 

Without Fence 

With Fence 

0 

0.589571 

0.59945 

0.091571 

0.0962 

3 

0.601254 

0.61923 

0.104455 

0.1069 

6 

0.61298 

0.62475 

0.1219 

0.1325 

9 

0.6313 

0.65841 

0.1228 

0.1395 

18 

0.72343 

0.73232 

0.128556 

0.1435 

19 


0.7593 


0.1484 



Figure 19: Variations of the Coefficient of 
Lift and Angle of Attack 
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0.16 
0.14 
0.12 
0.1 
rS 0-08 
0,06 
0.04 
0.02 
0 


Figure 20: Variations of the Coefficient of 
Drag and Angle of Attack 

6. CONCLUSIONS 

This project has successfully validated the increased effectiveness of delta wing by the use of wing fences. 
Through the CFD environment, the C L , C D and the angle of attack values for delta wing with and without wing fences were 
generated and a comparative study was conducted. This unrevealed the increased lift characteristics of the delta wing with 
fences. The average increase in C L is 13.03% with a small increment C D in say 9%. The stall angle has increased by 1.535°. 

The modification in aerodynamic characteristics of the wing leads to increase in the lift which tends to increase 
the aerodynamic load over the wing thereby leading to structural deformation. In order to ensure the modification done 
produces negligible structural deformation, the aerodynamic load only is taken into consideration and the Fluid-Structural 
Interaction (FSI) analysis is done, from which a positive result has been obtained that the structure deforms by a very 
negligible value say 5-8% than that of the actual deformation. 

In the real-time analysis a prototype was manufactured using Synthetic Laser Sintering (SLS) process using 
polyamide as material and was tested in a smoke tunnel, it was observed that the flow separation for delta wing with fences 
occurs only at a higher angle of attack than that without fences. It was also found that the delta wing with fences possesses 
higher stall angle than that without fences. 
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